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Abstract 

Using tools from the theory of random fields with stationary increments, 
we introduce a new class of processes which can be used as a model for the 
noise perturbing an SPDE. This type of noise (called harmonizable) is not 
necessarily Gaussian, but it includes the spatially homogeneous Gaussian 
noise introduced in [S], and the fractional noise considered in [3]. We 
derive some general conditions for the existence of a random field solution 
of a linear SPDE with harmonizable noise, under some mild conditions 
imposed on the Green function of the differential operator which appears 
in this equation. This methodology is applied to the study of the heat and 
wave equations (possibly replacing the Laplacian by one of its fractional 
powers), extending in this manner the results of to the case H < 1/2. 
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1 Introduction 

The theory of random processes, or random fields, with stationary increments 
was developed in the 1950's, the landmark references being [T^] and [55] ■ These 
processes have many interesting properties, most of which can be derived from 
their spectral representation as an integral with respect to a complex random 
measure with orthogonal increments. In particular, various classes of self-similar 
processes with stationary increments, also called "i?-sssi processes" [H being 
the index of self-similarity), received special attention. These include the Her- 
mite processes, which arise as limits in the non-central limit theorems (see [llj ) 
and the generalized grey Brownian motions (see [24]). The concept of operator 
scaling random field (OSRF) has been recently introduced as a replacement for 
self-similarity in the case of anisotropic fields. An explicit form of the spectral 
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density of a Gaussian OSRF with stationary increments was obtained in [7. An 
extension to the multivariate case was studied in |22| . 

The fractional Brownian motion (fBm) is the most commonly known process 
with stationary increments. There are two generalizations of the fBm in higher 
dimensions: the fractional Brownian field (fBf), and the fractional Brownian 
sheet (fBs). The fBf has stationary increments, being the only isotropic Gaussian 
i?-sssi process (see [23], [6]). The fBs was introduced in [18], and was studied 
by many authors (see [5T], [1]). It does not have stationary increments, but can 
be studied using spectral analysis due to its representation as in integral with 
respect to a complex Gaussian measure. In the recent years, the fBs has been 
used increasingly as a model for the noise in stochastic analysis (see [TJ], [T5]). 

The goal of this article is to introduce the tools necessary for the study of 
stochastic partial differential equations (SPDEs) driven by a noise which has a 
harmonic structure similar to that of a random field with stationary increments, 
in both space and time. For this, an important step is to define a stochastic 
integral with respect to the noise. 

In dimension d = 1, the definition of a stochastic (or Wiener) integral with 
respect to a process with stationary increments is usually taken for granted in 
the literature. In the case of the fBm, this problem was studied thoroughly in 
[23] and [in], but these references did not explore the fact that the spaces of 
Wiener integrands may contain distributions. Implicitly used by several authors, 
this fact has been recently proved in [TB] (in the case of the fBm) and [T7] (in 
the case of an arbitrary process with stationary increments). 

By analogy with the case d — 1, one may think that a natural space of 
Wiener integrands with respect to a random field {^ijtgRrf with stationary 
increments is the Hilbert space T-l defined as the completion of the set S of 
linear combinations of indicator functions with respect to the inner product 
(l[o,t], l[o,s])w = E(XtXs). However, in dimensions d> 2, even in the presence 
of the strong structure induced by the stationary increments, a collection of 
random variables indexed by the points in M.'^ is not rich enough to produce a 
Hilbert space H of integrands suitable for the study of a general SPDE. 

Therefore, in the present article we assume that the space-time noise is a 
collection X = {X{Lp);(p g T>c{^'^^^)} of random variables (indexed by the 
class I?c(R'^^^) of complex- valued infinitely differentiable functions on R'^+-'^ 
with compact support), which can be represented as stochastic integrals with 
respect to a complex random measure Ai with orthogonal increments on IR'^+^. 

In particular, the covariance of the noise is given by 

E[X{^)X{^)]^ [ J-^(r,0-FV'(r,6n(dT,dC), 

where J^tp denotes the Fourier transform of (f, and H is a symmetric tempered 
measure on R''"'"^, which is related to the control measure of A^. This definition 
produces a rich enough space H, and induces a structure which is amenable 
to spectral analysis, being very similar to the structure of a random field with 
stationary increments on R'*"''^. Moreover, this framework includes the Gaussian 
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noise considered by other references (e.g. [T], [2], [3], [5]) [S], [IH])) without 
being restricted to the Gaussian case. In particular, it provides a unifying 
framework for the study of equations with fractional noise in time, which covers 
simultaneously the cases H > 1/2 and H < 1/2. 

An interesting case arises when 11 = x /z, where is a tempered measure 
on R'', and i'((ir) = T~^i/(dr) is the spectral measure of a process {Zt)teR with 
stationary increments and covariance function R. In this case, the noise can 
be identified with a collection {Xt{(p);t > 0,(p g T>c{S.'^)} of random variables 
which admit the representation mentioned above and have covariance: 

E[Xt{ip)X,{ij)] ^ Rit,s)J{ip,i;), 

where J(<p,V') = /r^ •^'/'(O-^V^IOmC'^C)- 

In the Gaussian case, the process {Xt{(p); t > 0,(p ^ T>{R'^)} was considered 
as a model for the noise driving an SPDE by many references. These references 
considered only the case when the Fourier transform of in iS'(M'') is a locally 
integrable function / (or a measure F), and hence J can be represented as: 

J((p,V)=/ / ^{xMy)f{x-y)dxdy, V^,Vg5(R''). (1) 

jRd JRd 

See e.g. [8], [9], [10] for the case when Z is a Brownian motion, and [2]j [3]j [1] for 
the case when Z is a fBm of index H > 1/2. In general, the Fourier transform 
of /i in 5'(R'^) may not be a locally integrable function (or a measure). This 
is the case for instance, when d = 1 and fi{d(,) ~ |Cl " with a G (—1,0). The 
fact that representation ^ is not needed in the present article leads to a more 
general class of noise processes, even in the Gaussian case. 

This article is organized as follows. In Section[2]we review some background 
material related to random fields with stationary increments. In Section |3l 
we introduce the "harmonizable" noise X (described above), and we identify 
several classes of deterministic integrands with respect to this noise. In Section 
31 we define the random-field solution of a linear SPDE Lu = X driven by the 
harmonizable noise X, and we show that, under some mild regularity conditions 
imposed on the fundamental solution G of the equation Lu — 0, the necessary 
and sufficient condition for the existence of this solution is: 

2 

n(dr, < oo, (2) 

where FG{s, ■) denotes the Fourier transform of the distribution G(s, •) G 5'(M''). 

Section [5] and Section [6] are dedicated, respectively, to the parabolic case 
L = ^ — C, and the hyperbohc case L — — C, where C = — (— A)^/^ for 
some /3 > 0, assuming that U = i' x fi (with measures v and n as above). In the 
parabolic case, we show that condition ([2]) is equivalent to: 

Jri Jr + [1 + W(^)J-^ 



e-'^'J'G{s,-)iOds 
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where = |^|'^. In the hyperbohc case, condition ^ is equivalent to: 

/ , ^ I T- — T77^'^{dT)u(d£,) < oo. (4) 

When is the fBni of index H E (0, 1), conditions ([3]) and (|3]) become 



^(d^) < oo, 



extending the resuhs of ^ to the case H < 1/2. 

The appendix collects some auxiliary results, which are invoked in the article. 
In Appendix A, we review the construction of the stochastic integral with respect 
to a complex random measure with orthogonal increments, and list some of its 
properties. Appendix B contains some elementary estimates related to the fBm 
case. Appendix C shows that some technical conditions (which are imposed for 
treating the hyperbolic equation) are verified when the measure v is the Fourier 
transform of the Riesz kernel or the Bessel kernel. 

We conclude the introduction with a few words about the notation. We let 
'Dci^'^) be the set of complex- valued infinitely differentiable functions on R'' 
with compact support. For any p > 0, we denote by L^{R'^,fi) the space of 
complex- valued functions (p on such that \(p\P is integrable with respect to 
the measure ^. When ^ is the Lebesgue measure, we simply write L^{R'^). We 
denote by L^(S^) the set of complex-valued random variables X , defined on a 
probability space (f], J",P), such that E\X\p < oo. We let 5c(M'^) be the set 
of complex- valued rapidly decreasing infinitely differentiable functions on R"^. 
Similar notations are used for spaces of real- valued elements, with the subscript 
C omitted. 

We denote by J-^piO — J^d e~^^'''^ Lp{x)dx the Fourier transform of a function 
ip £ L^{M.'^). (We use the same notation for the Fourier transform of functions 
on R,]R'' or R'^+i. Whenever there is a risk of confusion, the notation will be 
clearly specified.) We denote hy x ■ y = ■^i^i ^^'^ inner product in R'* and 

by \x\ = {x ■ x)^!"^ the Euclidean norm in R''. A bounded rectangle in R** is a 
set of the form (x, y] = {z e R**; xi < Zi< yi for alH = 1, . . . , d}. We denote by 
TZd the class of bounded rectangles in R"* and by Bbi^"^) the class of bounded 
Borel sets in R"*. 



2 Background 

In this section, we review the basic elements of the theory of random processes 
and random fields with stationary increments. We refer the reader to Sections 
23 and 25.2 of [29] for more details. 

Definition 2.1 A real-valued random field X = (Art)jgjjd, defined on a prob- 
ability space (J7,J-", P), has (wide-sense) stationary increments, if E{Xf) < oo 
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for any t e M'', and the increments of the process have homogenous first and 
second order moments, i.e. E{Xt+s — Xg) = m{t) for any t,s £ M"*, and 

E\Xt+s - X,\^ = V{t) for any t,s e M''. (5) 

In particular, a stationary random field has stationary increments. 

In the present article, we work only with zero-mean processes, i.e. we assume 
that m(<) = for aU t e R'^. From ® , it follows that for any s,t,ue R'^, 

E[iXu+t - x^){Xu+s - = hy{t) + y(s) - Y{t - s)\ =-. R{t, s). 

(To see this, write {X^+t - X^){X.^+, - = ^{\X^+t - + |^«+. - 
X„p — jXu+t — X„+sp).) On the other hand, the covariance between arbitrary 
increments Xu+t — X„ and X„+s — Xy depends on t, s and u~ v: 

E[{Xu+t ~ Xu){Xy+s - Xy)] = R{t, s + v-u)- Kit, v-u). 

This follows by writing Xy+s - Xy ^ {Xy+s - Xy) - {Xy - Xy). 

The function R is called the covariance function of X and plays an important 
role in the analysis of the process X, via the spectral theory. More precisely, 
if R is continuous, then the process X and its covariance function admit the 
following spectral representations: (see (4.351) and (4.345) of _29j) 

Xt = Xa+t-Y+ f (e-*"* - l)7W(dT) (6) 
JrA{o} 

R{t,s) = T.t-s+ j {e'"-' -l){e"-' ~l)]l{dT). 

JR''\{0} 

The elements which appear in these representations are the following: 

• M denotes a collection {JV[[A)]A £ T?.^} of zero- mean random variables 
in Lp(r2), indexed by the class T?.^ of rectangles in R'^ which do not include 
0, which satisfy JJie following_conditions: 

(i) (additivity) M{AVJB) = M{A)+M{B) a.s. for disjoint sets A,B e 7^2; 
(a) (orthogonality) E[M{A)M{B)] = for any disjoint sets A,Be TZ^; 
(Hi) (symmetry) Ai{—A) — Ai{A) a.s. for any A e TZ'^; 
(iv) (control measure) E\M{A)\'^ = J1{A) for any A e TZ'^. 

• is a symmetric measure on M'^XjO} which satisfies the condition: 

/ T^^^idr) < oo. (7) 

The measure Jl is called the spectral measure of X. 

• y is a real- valued centered d-dimensional random vector with covari ance 
matrix S, which is orthogonal to (A) foraU A e 7^2, i-C- E[YM{A)] = 
for aU A e 7^°. We may define Mi{0}) = Y and J1{{0}) = E\Y\^. 
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Condition ([7]) ensures that the function (/?(t) = e~''^'* — 1 is in Lf.(R'',pt), 
and hence the integral of ([6]) is weh-defined. 

In the present article, we assume that Xq = and Y — 0. We consider an 
object Ai as above, except that its indexing class is extended to all rectangles 
(which may include 0), and the measure Jl is not required to satisfy ([7]). More 
precisely, we introduce the following definition. 

Definition 2.2 Let Jl he a symmetric measure on R''. A collection A4 = 
{A4{A); A E TZd} of zero-mean random variables in L'^{fl) which satisfy condi- 
tions (i)-(iv) above for any sets A,B E TZd is called a complex random measure 
on with orthogonal increments and control measure Jl. 

The stochastic integral with respect to A4 is defined as the mean-square limit 
of a sequence of Riemann-Stieltjes integrals. More precisely, for any function 
(fi (z L'^i^'^jjl) which is continuous on compact sets, we define the following 
C- valued random variable: 



where (^j)i<j<fc„ C TZd is a partition of [— -/V„, Nn] with Nn = {n, . . . ,n) £ R'*, 
\Aj\ is the length of the longest edge of Aj, and tj G Aj is arbitrary. For the 
sake of completeness, the details of this construction are given in Appendix A. 

We now examine the case d = 1. Let Z = {Zt)t^R be a real-valued centered 
random process with stationary increments and spectral representation ^ with 
Y = and spectral measure i>. We assume that Zq = 0. 

Let i/{dT) — T'^v{dT). Note that v \s & symmetric measure on R which 
satisfies the condition: 



The measure v is sometimes also called the spectral measure of Z. 
Since J"l[o.f](T) = er^'^ds = (e^*^* - l)/(-iT), it follows that: 

R{t,s) - E{ZtZ,) = f Tl[o.t]iT)T^o,s]irHdT). (9) 

Note that relation ^ does not hold in dimensions d> 2. 

Let S be the set of elementary functions on R, i.e. real linear combinations 
of indicator functions 1[q G R. We endow S with the inner product: 

(l[o,t], 1[o,s])a = Rit,s). 

Let A be the completion of S with respect to (•, •)a- The map l[o,t] >— > Zt is 
an isometry between S and L^(r2), which is extended to A. This extension is 
denoted by (^9 i— >■ Z{(p) = J^(p{t)dZt. The space A contains distributions. More 
precisely, if 93 G 5'(R) is such that is a function and Jj^ \J^(p{r)\'^i'{dT) < 00, 
then ip e A (see Theorem 3.2 of [17]). 





fc, 



lim { lim 

n->-oo maXj|Aj|->0 




(8) 
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Remark 2.3 By Theorem 2.3 of [T7] and Proposition 2.5 of [17], A coincides 
with the completion of I?(]R) with respect to the inner product: 



In general, the Fourier transform in iS'(M) of the tempered measure v is 
a distribution (denoted by J-v), which coincides with (\/2)V'\ V" being the 
second order distributional derivative of the continuous variance function V 
(see Section 2.2 of [IT]). In some cases, Fv is a locally integrable function p. 
By the Fourier inversion theorem on 5(M), this is equivalent to saying that v is 
the Fourier transform of (27r)^^/3 in 5'(]R). In this case, (see Lemma 5.6 of [ID]) 

R{t,s)= [ [ p{u-v)dudv. (10) 



^0 



Relation (|10p is useful in applications, but is not needed in the present article. 
When V is finite, p{t) = e~^'^^v(dT), and a process (Zt)tgK with stationary 

increments and covariance R given by pI7|) can be represented as Zt = /p Ygds, 
where (It)tgR is a zero-mean stationary process with EiYoYt) = p{t). 

Example 2.4 Let i^(dT) — \t\^^' for some 7 G (—1, 1). Then v is the Fourier 
transform of a locally integrable function if and only if 7 e (0, 1) (see Lemma 
3.1 of [IB])- This function is the Riesz kernel p{t) = c-y|f|~^^~'''' (see Lemma 1, 
Chapter V of [27j). In the Gaussian case, the parametrization 7 = 2H — 1 with 
H £ (0, 1) corresponds to the case when Z is a fBm of index H. 

Example 2.5 Let iy{dT) (1 + I'rp)^'''/^ for some 7 > —1. When 7 > 0, is 
the Fourier transform of the Bessel kernel (see Proposition 2, Chapter V of [77] ) 

/■oo 

p{t) = c, / u;(^-i)/2-ie-'"e-*'/(4-)d^z;. 



The measure is finite if and only if 7 > 1. If 7 = 2, then p(t) = ci exp(— C2|i|) 
for some constants ci,C2 > (see Exercise 2.3.4, Chapter 8 of [E]), and 
a process {Zt)tem. with stationary increments and spectral measure vidr) = 
T~^(l + T^)^^(ir can be represented as Zt — J^Ygds, where (It)tgR is a zero- 
mean stationary process with E(YoYt) — p(t). In the Gaussian case, (Yt)teR is 
a zero-mean Ornstein-Uhlenbeck process. 

3 The harmonizable noise 

Let n be a (symmetric) tempered measure on E'*+^, i.e. 

/ f -, , 2 , 1^12 ) n(dT,dO<oo for some fc>0. (11) 
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Let TW = {A4{A); A E TZd+i} be a complex random measure on R'^+-'^, with 
orthogonal increments and control measure 11 defined by: 



n(dr, dO = -YPT- — 72 n(c^T, for r e M, ^ e 



For every S 5c (R'^"''^), we define a C- valued random variable by: 



d 



Xi^)= f T^iT,OHr)fli-iQMidT,dO. (12) 

To see that X{(p) is well-defined in L'^{il), we use the fact that Ttp G 
5c(IR'*+i). Hence sup^gK,«6R'i(l + + l^l^TlJ^V^ir, < oo for any ti > 0. 
Taking n = k and using ([lT|) . it follows that 

i.e. the function i/'(r,0 = -^<^(t, nj=i(-«0) belongs to L^. (Rd+i^ fj). 

Since ■0 is also continuous on compacts, X{(p) is well-defined (see Appendix E)) . 

If is real- valued, then X{(p) is real- valued, due to the symmetry of A4 (see 
relation ([75]) . Appendix A). Using ([7^ . we see that for any ip,^p G Sc 



E[X{ip)Xi,p)] = / T^{t, 0-FV'(t, on(dT, =: I(^, V)- 

The functional I is bilinear, symmetric, and non-negative definite. To see 
this, note that for any m > 1, ai, . . . , a„i G C and . . . , t/^m e 5c (R''^"'^), 



m „ 



i=l 



2 

n(dr, dO > 0. 



We identify any functions (p,ip € 5c(R'*^"'^) for which X{ip — tpjip — ^p) = 0, 
so that I is an inner product. 

Remark 3.1 If the Fourier transform of 11 in 5'(M'*+^) is a locally integrable 
function F : R'^+i -)■ [0,cx)), then 

X{ip,^p) — / / ip{t, x)^p{s,y)F{t — s, X — y)dxdydtds, 

for any ip,ip G iS(R''^^). This representation is not needed in the present work. 

In what follows, we work with functions (p e I?c(K'*^^). More precisely, we 
introduce the following definition. 
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Definition 3.2 The process {X{(p); ip G 'Dc{M.'^+^)} defined by ^ is called a 
harmonizable noise with driving measure M and covariance measure 11. 

We endow I?c(R''^^) with the inner product (•, defined by: 

and we let H be the completion of Vci^'^^^) with respect to this inner product. 

The map (p i~> X{lp) is an isometry from 'Dc{M.'^^^) to which is 

extended to T-l. The Hilbert space H contain distributions. By abuse of notation, 
we write 

X{ip)^ / / (p{t,x)X{dt,dx) 

even if 1^9 is a distribution. This stochastic integral is well-defined only ii ip £ H. 

We will need an alternative representation for the inner product in H. For 
any (pi,(p2 G 'DcO^'^^'^), if we denote by 

the Fourier transform of the function (pi{t, •), and by ^(j)^^^ the Fourier transform 
of the function 1 1-^ (0 for i — 1,2, then 

{^i,P2)n= f Tcj)fHT)T^f{T)n{dT,dO^:{m,V2)o- (13) 

The following result is a modified version of Lemma 3.7 of [2]. 

Lemma 3.3 Let ^ be a tempered measure on M . For every A G Bb(W^), there 
exists a sequence {ipn)n>i C I?c(]^'^) such that 

Proof: As in the proof of Lemma 3.7 of [2], let r]n{x) = n'^ri{nx), where rj G 
Vci^'^) is such that rj > and / r]{x)dx = 1. Setting (j)^ = 1a* r]n, we note 
that 0„ e 'Dci^'^), 0,1 — ?► 1a uniformly, and supp (j)n C K for all n, where K is 
a compact set in R''. By the dominated convergence theorem, 

Jr-^ 

For the second part of the proof, we argue as on p. 613 of [17]. Since the 
Fourier transform is a homeomorphism from onto itself, and I?c(K'') is 

dense in 5(R''), J"(X>c(R'')) is dense in 5c(R''). Since /i is a tempered measure, 
the convergence in iSc(IR'') implies the convergence in L'f.{W^,fi). Hence, there 
exists a sequence (V'n)n>i C 2?c(R'') such that: 

/ IMO-J'MOl'KdO^o. 

Jr'' 
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The conclusion follows. □ 

Next, we give a criterion for checking that ip Cz H, which is a generalization 
of Theorem 2.1 of [3]. 

Theorem 3.4 Let M. 3 t i-^ (p{t, ■) e S'{M.'^) be a deterministic function such 

that J-(p{t, •) is a function for all t €E R. Suppose that: 

(i) for all £, € M'*, the function 1 1->- Tip{t, •)(^) =: (t>^{t) belongs to 

(a) the function (t,^) i— >■ J'fit, ■){£,) is measurable on R x R'^. 

For any ^ G R'', we denote by the Fourier transform of (j)^. If 



l^llo- / |J-0e(r)|2n(dr,dO <oo, 



(14) 



then ip G H and W^fiW^ = H^'Ho- -^'^ particular, the stochastic integral of ip with 
respect to X is well-defined and is given by: 



r. d 

/ T^^{T)i-tT)l[{-t^,)MidT,dO. (15) 



Proof: The argument is a simplified version of the proof of Theorem 2.1 of [3]. 
Due to (IT51) and the definition oi H, we have to show that for any e > there 
exists l = le& Vc{W^"^^) such that 

\W^l\\l^ I \F^^{T)-Fij^{T)\'n{dT,d£,)<e\ (16) 

where F4>^ denotes the Fourier transform of the function t ^ ip(^{t) — J- l{t, ■){£,). 
We denote a(T, ^) — J^0^(t). It will be shown later that: 

the function (t, ^) i-^ C) is measurable. (17) 

Due to (HH) and ([HI), it follows that a e L^(R''+\ 11). By applying Theorem 
19.2 of [S] to the real and imaginary parts of the function a, we infer that there 
exists a simple function h on R'^+^ such that 



/ |a(r,0-M^,OI'n(dr,de)<^. 



(18) 



Here, a simple function is a complex linear combination of indicator functions 
of the form 1b, where i? is a Borel subset of R''^^ which is not necessarily 
bounded. From the proof of Theorem 19.2 of fS^, we see that the Borel sets B 
which appear in the definition of the function h satisfy Tl{B) < oo. The measure 
n is locally finite, and hence cr-finite. By applying Theorem 11.4.(ii) of 5^, each 
of these Borel sets can be approximated by a finite union of bounded rectangles 
in R'^"'""'^. Therefore, we may assume that each set B is bounded. 

By Lemma [3751 there exists a function I E I?c(]R.''^^) such that 

/ \lB{T,0-n{r,0?UidT,dO < ^. (19) 
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Relation (dH) follows from ^ and ((H]), since J'Z(r,0 = J^ipd'^)- 
We now prove ([T7)) . Since (t, i->- is measurable, by applying Theorem 
13.5 of [5] to the real and imaginary part of (j)^{t) we conclude that there exists 
a sequence (/„)„ of complex-valued simple functions such that ln{t,(,) — </'j(t) 
and \l„{t,^)\ < |<^4(i)| for any n > 1. We let a;„(r,0 = J'?„(-,e)M- By the 
dominated convergence theorem, 



oo 

ITt ( 



Jo 

If /(t, ^) = l(c.d](Ol^(C) is an elementary function, then ai{T, ^) = J-l{-, ^){t) = 
•^1(c,(/](''')1a(C) is clearly measurable. Since ai„(r, ^) is measurable for any n, it 
follows that a(r, ^) is measurable. 

To prove (jlSp . we denote by y the random variable on the right hand side of 
P5|) . which is well-defined due to and hypothesis fij. By definition, = 
lim„_s.oo -'^(Vti) in L(.(r2), where ((^n) C 'DcQS.'^'^^) is such that ||(p„ — (^||o — > 0. 

Using ((HI) and ([71|) (Appendix A), we have: 



E\Xiipn)-Y\'' = E 



|^^„(r,e) - J-0e(r)pn(dT,dC) = llv'n - ^llo ^ 0, 

-1 

i.e. Y = lim„_s.oo X{(pn) in L'^{il). Hence, Y — X{lp) a.s. □ 

The next result gives a criterion for checking if the indicator function of a 
bounded Borel set A (in particular, a bounded rectangle) is in H, and gives a 
representation for the stochastic integral of 1a- 

Theorem 3.5 (i) If A e Bb{R'^+^) is such that 

f |j-u(T,oi'n(dT,dC) <(X3, (20) 

then 1a (zH and the stochastic integral of 1a with respect to X is given by: 

X{A) X{1a) = / J-U(r,0(-*T) \\{-ii,)M{dT,dO 
(ii) Suppose that the measure H satisfies the following condition: 

Then I120\) holds for any finite union A of bounded rectangles in R'^^^. In par- 
ticular, the stochastic integral o/ l(o,t]x(o,a;] ™th respect to X is given by: 

P d 

X(i,x):=X(l(o,t]x(o,.])= / (e-^^*-l)n(e^'^^"^-l)M«fT,rfO- (22) 
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Proof: (i) Let ~ 1a * (pn be the mollification of l^i, where (j)n{t,x) = 
n'^+^(t){nt,nx), and </> e 'D{W^+'^) is such that |!(?!)||Li(Rd+i) = 1. Then (p„ G 
V{W^+^), Fipn - :FIa = -FlA(-7^0n - 1) ^ and - 1| < 2. By the 

dominated convergence theorem, 

/ |^^„(t, - ^u(t, OI'n(dT, 0. 

This proves that 1^ S The second statement follows by an approximation 
argument, as in the last part of the proof of Theorem 13.41 

(ii) By linearity, it is enough to consider the case when A = (0, t] x (0, x] with 
t G G M'^. Note that ^ is equivalent to A (t^ + |^|2)]n(dr, < 

00. On the region R = {(r,C) G M^'+^r^ + j^p < 1}, we use the fact that 
< 1^1, where |^| is the Lebesgue measure of A, and hence 

/ \TlA{T,0\MdT,0 < / U{dT,0 < oo. 

JR Jr 

On the region R'^, we use the fact that 



l-^U(r,OI - 



-IT 



n 



< 2 



d+1 



1 



and hence, 



□ 



|J"lA(T,OPn(dT,dO < 2^+1 / n(dT,dO < oo 



When condition (I^TI) holds, there is a one-to-one correspondence between 
the class of random fields a;)}(j ^.jgjjd+i defined by and the class of 

random fields with stationary increments and spectral representation: 



X*{t,x) 



(23) 



This correspondence is illustrated by the next two examples. 



Example 3.6 The fractional Brownian sheet (fBs) with indices H, Hi, . . . , G 



(0, 1) is a zero-mean Gaussian random field {X{t, x)}tt 



-1 with covariance 



E[Xit, x)Xis, y)] = Rnit, s) J] Rh, {x,,yj), 



where Rnit, s) — {\t 



2H 



\2H 



t — s| )/2 is the covariance of the fBm of 



index H G (0, 1). Since 

RH{t, s)^CH f [e-"' - l){e"^ - 1)\t\-^^"+'Ut, 
Jr 
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with ch = T{2H + 1) sin(7riJ)/(27r), the fBs has the representation (P^ . where 
M is Gaussian and n(dr, = CH|r|-(2ff+i) lf.^^{cH^ 101"^^^^+^^)- The mea- 
sure n satisfies (PT|) if and only if _ff + X]j=i < 1- this case, there is an 
associated random field with stationary increments defined by which is 

also called a "fractional Brownian sheet" by some authors (see T)- 

Example 3.7 A fractional Brownian field (fBf) on with index H € (0, 1) 
is a zero-mean Gaussian random field a;)}(f ,j.)gRd+i with covariance 

E[X*it,x)X*is,y)] = i[(t2 + \x\T + is' + \y\Y {\t ^sj^ + \x - y\T]- 



This random field has stationary increments and spectral representation (|23p . 
where M is Gaussian and U{dT,d^) = (t^ -|- |^|2)-(2H+d+i)/2 rpj^g measure fi 
satisfies condition ((2T|) for any e (0, 1). To this process, one may associate a 
random field {X{t,x)} defined by 

An interesting case arises when 

n = iyx ^, (24) 

where is a symmetric measure on M satisfying ([5]), and /i is a symmetric 
tempered measure on M.'^. Using the same mollification argument as in the 
proof of Theorem [X51(i). it follows that l[Q.t](p G "H for uny t£R,(p e VcO^'^)- 
To see this, note that 

\HMo.t]V>)ir,0\'mr,dO= f "'~^^\ idT) [ \^ipiO\'fi{dO<oo. 

For any t,s E S. and E 'Dc{R'^), we have: 

(l[0.t]<y5, l[O,s]'0)« = R{t,s)J{(p,1p), 

where R is given by ([9]), being the covariance of a process Z = (Zf)tgR with 
stationary increments and spectral measure vidr) — T^'^v^dr), and 



In this case, T-l is the completion of £ with respect to the inner product 
(■)■)«) where £ is the set of complex linear combinations of functions of the 
form l[o,t]<i5, with < e M and (p E 'Dc{R'^). To prove this, we use that fact 
that a function in 2?c(K''+^) 

can be approximated by a function of the form 
(p{t,x) = ipi{t)ip2ix) with ipi E X'c(K),<^2 e VciM."^), and a function in VdM.) 
can approximated in the norm || • ||a by a complex elementary function on M 
(see Remark [23|) . 

The stochastic integral of l[o,t]'f with respect to X is given by: 

Xt{^) := X{l[o,t]^) = / (e— ' - 1)J^^{0 Y[{-ti,)M{dT,dO. 
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4 Equations with harmonizable noise 

As in the previous section, let {X{ip);(p G 2?c(IR.''+^)} be a harmonizable noise 
with driving measure Ai and covariance measure H. We consider the equation: 

Lu{t,x) ^ X{t,x) t>0,xeR'^ (25) 

with some deterministic initial conditions, where L is a second order partial 
differential operator. Let w be the solution of the equation Lu = with the 
same initial conditions as (1^51) . and G be the fundamental solution of Lu — 0. 
The following definition introduces the concept of solution. 

Definition 4.1 The process {u{t, x);t > 0,x E M''} defined by 

u{t,x) =w{t,x) + [ [ G{t- s,x-y)X{ds,dy) (26) 
Jo Jw 

is called a random field solution of (|25p , provided that the stochastic integral on 
the right-hand side of is well-defined, i.e. for any t > and x e K'^ 

9tx := l[o,t]G(i - •, a; - •) e "H. 

The following result gives some general conditions which ensure the existence 
of a random field solution, as a direct consequence of Theorem 13.41 

Theorem 4.2 Assume that G{t,-) e 5'(R'^), such that TG{t,-) is a function 
for any t > 0, (t,^) TG{t, •)(^) =: H^{t) is measurable on R+ x R'', and 

[ \TG{s, ■){^)\ds < oo, for anyt>0,^e'R'^. 
Jo 

(i) Equation V25\] has a random field solution {u(f, x)}(j if and only if 

It := / \Fa.tH^{T)\^W{dT,dC) < oo for all t > 0, (27) 

where J^o,tH^{T) = e^'^^''H^{s)ds, r e M. In this case, E(u{t,x)) ~ w{t,x), 
Var(it(i, a;)) = It and the solution is given by: 

u{t,x) = wit,x) + / e-^^*-'«•"J■o,ti^^(r)(-^T) T\{-tCj)M{dT,dO. 

(a) Assume zero initial conditions. Then, for any t, s > and x,y E M.'^ , 

E[u{t,x)u{s,y)]= f e-'"(*-^)e^'«-("-^)j-o,t775(T)J-o,«i/5(T)n(dT,dO- 
In particular, for any t > fixed, 

E[u{t,x)u{t,y)]^ f e-'i<^-y^fitidO. 

where ^t{) =nt(Rx-) andIlt{dT, d(,) = |Jx),t^^c(''')Pn(dr, d^), i.e. {u{t,x)}^(zRd 
is a zero-mean stationary random field on with spectral measure ^t ■ 
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Proof: (i) We apply Theoreni l3.4l to — gtx- Note that the function 



-.^TgU^-m = l[o,t](s)e-'«-"^G(t-s,.)(0 = l[o.t] (s)e-*«-i7^ (t - s) 



satisfies conditions (i) and (ii) of Theorem 13.41 

Let J-cj)^ be the Fourier transform of the function s i— > (f>^{s). Then, 

J-(/)^(t) = J e-*"(/)^(s)ds = e^*«-^^ e""' H^{t - s)ds 

= e-*«-^e-*^* / e"'l^(s)ds = e-^^'^e"^^* JoT^TW- (28) 

"'0 

By definition u(t,x) = w{t,x) + X{gtx)- Hence, E(u(t,x)) — w{t,x) and 

Y'^r{u{t,x)) = E\X{gtx)? = \\gtxfu=If 

(ii) Using the isometry property of the stochastic integral i— >■ X[ip) and 
the fact that WpWu = where || • ||o is defined by relation ([14]), we obtain: 

E[u{t,x)u{s,y)] = E[X{gtx)X{gsy)] = {gtx,gsy)H = {9tx,gsy)o 



F^^{T)F4,^{T)T\{dT,di), 

-1 

where (/'{(u) = J-gtx(u, ■)(^) and tp^iu) — J-gsy{u, ■){^). The conclusion follows 
using (gll). □ 

5 A parabolic equation 

In this section, we assume that ([M]) holds. We consider the equation: 

-^{t,x) + {~A)^^^uit,x) = X{t,x), t>0,xeR'^, (29) 
for some /3 > 0. In this case, 

H^it) ^ TG{t,-) ^ exp{-WiO}, 

where 5'(^) — c^lCI^ and > is a constant. By Theorem 14.21 it suffices to 
check that p7)) holds. For this, we write /( = Jj^^ ^tiOf^i^O: where 

NtiO ^ ljJ'o^tH^{r)\'HdT). (30) 

We proceed to the calculation of Fo^tH(^{T). For this, note that if ip{x) = 
e-°^ with a e M, then 



"'0 



1 - e-'^°-+"^^ _ 1 - e^'^* cos(rt) + ie^"* sin(Ti) 
a + ir a + ir 

and 

\^oMr)\^ = ^^{sin2(rt) + [e-'^* - cos(rt)]2}. (31) 
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Applying (PT|) with a = ^'(^), we obtain: 

l-^o,ti^e(r)P = ;^^^^{sin^(rt) + [e-**(«) - cos(rt)]^}. (32) 
Using ([50)1 and ([5^ . we obtain: 

To identify the necessary and sufficient condition for the existence of a solu- 
tion, it suffices to obtain some suitable estimates for Nt{^). This goal is achieved 
in the next theorem. 



Theorem 5.1 For any t > and ^ e M'', we h 



ave: 



Jr^ r2 + [1 + ^(O]^ ^ ' - - * t2 + [1 + ^'(e)]2 

where and Cj^'' are some positive constants depending on t. Consequently, 
equation i29\) has a random field solution if and only if (0j holds. 

Remark 5.2 In the case of Example 12.41 and Example 12.51 with 7 e (—1,1), 
condition ^ becomes: (use Lemma [B. 11 Appendix B with a = 1 + 5'(^)) 

I \ 1+7 



fi{d£,) < 00. (33) 

When applied to Example 12.41 with the parametrization 7 = 2H —1,Hg (0,1), 
Theorem 15. II becomes an extension of Theorem 3.4 of [I] to the case H < 1/2. 

Proof of Theorem 15. It (a) We first show the upper bound. We claim that: 

NtiO < 2iv:max(t^ 5), for aU ( e R"^. (34) 

To prove (|34p . we split the integral on the right-hand side of (15(7)) into the 
regions |t| < 1 and |r| > 1. We denote the two integrals by Nj:^\s^) and N^'^\(_). 
For Nf ^\^), we use the fact that: 



e 



< / e-^*(«)ds<i, 



and hence 

N^'\0 = I \T,MT)\''dj.{T) <t^ I ^{dr) <t' [ -^iy{dr). 

J\t\<1 J\t\<1 •'\r\<l -'-+''' 

(35) 

For {^), we use the fact that: 



snr 



'■{Tt) + [e-**(«) - cos(rt)]2 < 5, (36) 
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and hence 

i|T|>l ^ "'|t|>1 ^ "'|r|>l ^ + 

Relation ^ follows from §5^ and dSZ]). On the other hand, by 

iV«(0<5^ ^^^|^^^^^ ^(dr), foraUeeM'^. (38) 
('jj Assume ^'(^) < 1. Then + [l + *(^)]2 < +4 < 4(r2 + 1), and hence 
^ i^(dr) > i / -^jy(dr) = ^K, if ^(O < 1. (39) 



/r t2 + [1 + Vl/(^)]2 V ^-47^1 + ^2 V ' 4 

From ([M]) and dH]), we infer that: 

1 



NtiO < 8max(^^5) -^-^^-^^^.(dr) if ^{0 < 1- (40) 

(ii) Assume ^'(^) > 1. Then 4[t2 + *(^)]2 > 4t2 + [1 + *(0]^ > + [1 + 
\I'(^)]^, and hence 

From ([38]) and (gl]), we infer that: 

m)< 20 l^^^^^^^Hdrl iinO>l- (42) 

The upper bound follows from (gO]) and (|42|), with C^"^^ = 8max(i2,5). 

(b) We now show the lower bound. Let pt > be a constant depending on 
t, which will be specified below. 

(i) Assume 5'(^) < pt- We use the fact that sinx is decreasing on [§, tt]. Let 
|<c<c?<7rbe fixed. We have: 



C<|Tt|<d 



^2+^(^)2- 



> sin^ d I r2~T~a777T2 '^(^''') 

> ,^2 / Krfr) =: Ait). (43) 



d^/t^ + Pi Jc/t<\T\<d/t 

On the other hand 
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From (132]) and (gH, we obtain: 

('iij Assume ^(^) > p*. For any x > 0, let 

^t(x) = sin2(rt) + [e-*=" - cos(r<)]2 = 1 - pe'*^ cos(rt) - e-^*^]. 

Note that Jt,T{x) — >■ 1 as x — >■ oo, uniformly in t G M. Hence, there exists pt > 
such that 

fr,t{x)>-, for all X > and T G R. (46) 

(More precisely, |2e-*"= cos(Ti) - e^^t^] < 2e-*^ + e^^te < 2c + c^, if e"*^ < c. 
Choose c G (0, 1) such that 2c + < 1/2, e.g. c = 0.1. Then /t,r(x) > 1/2 for 
any x> {— lnc)/t =: pt-) Using P5|) . we infer that: 

^ ^ L ^ ^ I .- + [iU(e)F -(^-)' " 

(47) 

The lower bound follows from (gS]) and (gT]), with C^*^^ = min i|. □ 

Remark 5.3 The results of this section are also valid for the equation 

^{t,x) - Cu{t,x) = Xit,x), t>0,xe R"^, 

where £ is the L^(M'')-generator of a d-dimensional Levy process {Xt)t>o with 
(real-valued) characteristic exponent 5'(^) (as in 13!). Assuming that Xt has 
density pt, we see that G{t, x) = pt{—x) and 

i?4(<) = FG{t, -m = / e'i-^pt{x)dx = £;(e*«-^') = exp{-i*(0}. 

6 A hyperbolic equation 

As in Section [5l we assume that (l24l) holds. We consider the equation: 

-^it,x) + {-Af/^u{t,x) = X{t,x), t>0,xeR'^, (48) 
for some /3 > 0. In this case. 



where ^'(0 = (see (16) of [10]). 
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Following the same steps as in the parabolic case, it suffices to find appro- 
priate bounds for It = J^^d Nt{£,)fi{d^), where Nt{£,) is given by (|30)) . For the 
calculation of J^o,t-ffc (t), we use the change of variable r = s-y/w(^: 

ToMt) = -i= Te— ^sin(sv/*(|))ds (49) 
1 



1 



tV*(«) , 

g-»rr/V*(«) sin rdr 



J'o.T sin , where T = t V*(0- 



An elementary calculation shows that: (see the proof of Lemma B.l of [5]) 
|J-o.Tsin(r)p = ^-^^[f^{r)+gUr)]. (50) 

where /t(''") — sin(TT) — rsinT and griT) = cos(rr) — cosT. Hence, 

[sin(rt) - sin(tv/¥(|))]2 + [cos(ri) - cos(iVW))]' 

\^oMr)\'- = (^2 _ viy(^))2 • 

(51) 

Using (|30|) and (|5T1) . we obtain: 



= /. — 

We assume that 

i/(dr) = ?7(|r|)dT, 
where the function rj satisfies the following condition: 

(C) for any A > there exists C\ > such that r]{\T) < Cxrjir) Vr > 0. 

In addition, we assume that rj satisfies either (CI) or (C2) below: 

(CI) rj is non-increasing on (0, oo), and for any K > there exists Dk > 

such that Jq rj{T)dT < DKari{a) for any a> K 
(C2) rj is non-decreasing on (0, oo), and for any K > there exists Dk > 

such that T~^ri{T)dT < D/f a^^ry(a) for any a > K. 

Remark 6.1 In Examples 12.41 and 12.51 condition (C) holds for any 7 G (— 1, 1), 
respectively for any 7 > — 1. In these two examples, condition (CI) holds if 
7 G [0, 1), whereas condition (C2) holds if 7 € (—1,0) (see Appendix C). 

Under these conditions, we obtain the following result. 
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Theorem 6.2 Suppose that-q satisfies (C). In addition, suppose that-q satisfies 
either (CI) or (C2). Then for any t > and ( S 

where D^^'^ and -D^^'' are positive constants depending on t, and 



N{C) - , / yjdT). 

Consequently, equation has a random field solution if and only if 

-v{dT)ii{d£,) < oo. (52) 



VT+m {\r\ + v/r+W))^ 

Remark 6.3 Note that condition ((52|) is equivalent to (|4]). This fohows using 
the fact that r^+b^ < (|t| + b)^ < 2{t'^ + b^) with b = ^l + *(0. 

The proof of the theorem foUows from the lemmas below. The first two 
lemmas treat the upper bound. 

Lemma 6.4 Assume that rj satisfies ( C) and is either non-increasing on (0, oo), 
or non- decreasing on (0, oo). Then for any t > and ^ G , 

NtiO < Ct^l= [ — i==-z.(dr), 



VmJ^i\r\ + vw)r 

where C is a positive constant. 



Proof: Using the notation a — we have: 



(r^ — 

We denote by N^^\^), Nf^\^) and Nj:^\^) the integrals over the regions 
Ri = {\t\ < a/2}, i?2 = {|r| > 3a/2}, respectively R3 = {a/2 < \t\ < 3a/2}. 
(i) We first treat nI^\£/) and n[^\s/). For this, we use the inequality: 

/|(t) + gl{T) < 2T{1 + \T\f for any r G R. (54) 

(To see this, note that |/t(t)| < 1 + |t| and |/t(t)| < 2T\t\, since | sina;| < |a;|. 
Hence, /2(r) < 2T|r|(l + |t|). Similarly, \gT{T)\ < 2 and |5t(t)| < T(l + |r|), 
since |1 - cosa;| < |a;|. Hence, g^ir) < 2T(1 + |t|).) 
Using ([54)) . we obtain: 

t\V 2t 



fL(-)+gl(-)<2ta(l+^-^) ^-{\r\+ar. (55) 
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Therefore, 



a Jr, {\t\ - ay 



for f = l,2. 



When r e both (|t| — a) ^ and (|t| + a) ^ behave as a ^, since 



< 



4 , 4 
< — and - — - < 



< 



a2 - (|r|-a)2 " 9a^ " (|r|+a)2 

When r e i?2, both (|t| — a)^^ and (|t| + a)^^ behave as t"^^ since 



1 



< 



1 



9 

< — and 



9 



< 



1 



T-' (|T|-aj^ 25t2 (|r| + a)^ 

Letting C(i) = 18 and C^^) = 50, we obtain 



< 



1 



ajR, (|t| +a)' 



:v{dT), for i — 1,2. 



(a) We now treat N^^\0- From (gH]) and we see that: 



(t2 - o2)2 



fta ( „ 



(56) 



Assume first that rj is non- increasing on (0,cxi). Then 



N^'\0 = ^ I \J^o.tSm{a ■){r)\M\r\)dT 



1 



< 



via/2) 



|J"o,tsin(a ■){T)\^dT 



via/2) , /•* . 2. ^ , ^ ^ ,^(« 
= ^ — 2n / sm (asjcts < znt ; 

a Jo « 



2 ■ 



(57) 



where we used Plancherel's theorem for the second equality above. On the other 
hand, when r G R3, (|r| + a)~2 behaves as a~^, since: 



< 



1 



< 



25a2 - (|t| +a)2 " 9a^ 
Using the fact that 77 is non-increasing, we obtain: 

^ I TT^TT^dr > - ■ ^vCia/2) 

IR3 



From condition (C), there exists a constant C1/3 > such that: 

?7(a/2) < Ci/3?7(3a/2). 
From dSD), dSHl) and dSH), it follows that: 



(58) 



(59) 
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When 77 is non-decreasing, the argument is similar. In this case, 
,(3),.^ / o_.^(3«/2) / V{al2) ^ 25 1 /• 1 



8 a J,^^ (|r| + ay 



□ 



Lemma 6.5 Assume that rj satisfies (C) and is either non-increasing on (0,oo), 
or non- decreasing on (0, 00). Then for any t > and ^ G Mf^ , 

Nt{0<Di''>^=l== [ — -L==-,.{dr), 



where D^^^ is a positive constant depending on t 



Proof: We denote a ^ b = ^1 + and 

^^^^ = / 7]~n~^^(^'^)' x>o. 

(i) Assume that a < 1. We first show that: 

NtiO < 3if max{^,8(t2 + 1)}. (60) 

Using ((53)) . we denote by Nf^\s^) and Nf'^\^) the integrals over the regions 
|r| < \/2, respectively |t| > \/2. Since | sinx| < a; for all a; > 0, we obtain: 

1,^ . , M /■* |sin(as)| , f' 
a Jo a, Jo 2 

and hence, by 

■/|r|<\/2 a 4 J|r|<V2 

r 1 

'|r|<V2 

When |t| > a/2, we have — > ^r^, since a^ <\ <\t'^. Using the fact 
that |/T(-r)| < 2T\t\ and |5t(t)| < 2, we obtain: 

fl (-) + .9?a (-) < 4i^«^4 + 4 = 4(tV2 + 1). 

Hence, 

i\ f 1 



^ ^ / (61) 



< 24(i2 + l) / ^^Krfr). (62) 
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Relation (|B(71) follows by taking the sum of (IFTj) and (p^ . 

On the other hand, since 1 < 6 < \/2, (|t| + 6)^ < (|t| + \/2)2 < 4(t2 + 1), 
and hence 

From dnO]) and we obtain: 

< 12v^max{^, 8(^2 + =: A^. (64) 

(ii) Assume that a > 1. In this case, a < b < ^/2a. By Lemma we have: 

1 /2 
NtiO < Ct-I{a) < Ct^I{a). (65) 
a h 

For any a; > 0, we write I{x) — Ii{x) + l2{x), where 
Note that 



h{a) < ^ / .{dr) < 1 / 

« J|r|<a J\t\ 



l\T\<a ^ ■'\r\<b 

Using condition (C), we have: 



h{a) < [ ^HdT)<[ J,v{\r\)dr ^ V2 [ ^ ^ ( ^ |r|) dr 
< ^/2Ci/^ / \y{dT) < 4V2Ci/^/2(6). 

J\T\>b ^ 



Therefore, 

J(a) < max{8,4\/2Ci/^}/(6). (66) 
From ((65)) and (166]) . we conclude that: 



Nt{i) < C%/2max{8,4V2C7i/^}ti/(5) =: ^i^- (67) 

The conclusion follows from and ([57)) with Dj^-* = maxjUj, □ 
We now treat the lower bound. 

Lemma 6.6 Assume that rj satisfies ( C), and is either non-increasing on (0, oo), 
or non- decreasing on (0, oo). Then for any t > and ^ € M.'^ , we have: 



Nt{0>C't^^^^^, tf ty^)>l, 



where C is a positive constant. 
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Proof: Let a = ^/'^iO- Assume first that 77 is non- increasing on (0,oo). 

Let e > be arbitrary (to be chosen later). Since the integrand of (|53|) is non- 
negative, Nt{^) is bounded below by the integral over the region {|r| < (l+e)a}. 
In this region, using the fact that 77 is non-increasing and condition (C), we have 



v{\t\) > J7((l +£)a) > C^/(i+^)rKa) =: c,7]{a). 



Hence, 
where 

h 

h 



\T\<{l+e) 



fta ( _ 



■9- [a 



dr =: c. 77(a) (/i- 12), 



(r2 - a^y 



fta I 



at 



dr 



\T\>{l+e)a (t^ -a^Y 



fta \ _ 



9t 



''id 



dr. 



Using Plancherel's theorem and the fact that 1 — (sina;)/x > 1/2 for any x > 2, 
it follows that: 



h = 



[ \J^o tsm(a AirM'^dr — —271 [ sm^(as)ds 
1 Jm Jo 

sin(2at) 



1 - 



2at 



IT t 

2^ 



On the other hand, using (j55p . we obtain that: 



2t 



a J\r\>(l+e)a (kl " af 

Choose e > such that f > p Hence, 



dr ~ — / -^dr = ^. 



TT 4 

2 ~ £ 



If 77 is non-decreasing on (0, 00), we bound Nt{£,) below by the integral over 
the region |t| < (1 — e)a, for some e G (0, 1). In this region, 

V{\r\) > r/((l - e)a) > C-^\,_,^via) 4?7(«)- 
As above, we obtain that iVt(^) > c^7/(a)(/i — I2), where 

1 



\T\<{l-e)a -a2)2 
1 



2 
9ta 



dr 



|r|<(i-£)Q (a - h\y 



TdT 



At r 1 



-dr = 4 
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Choose £ G (0, 1) such that f > 4(i - 1), i.e. !>£>(! + Hence, 
□ 

Lemma 6.7 Suppose that 77 satisfies either (CI) or ( C2). Then for any K > Q 
there exists a constant Mk > such that 

~ [ 7n~, — ^^{^'^) ^ Mk^^^^, for any a> K. 
a Jr [\t\ + ay 

Proof: If (CI) holds, then 

77(a) 



aJ\r\<a[n+ay Jo 



a 



l\T\<a 

1 /■ 1 /, ,^ , 277(a) /""^ 1 , „77(a) 

J\T\>a{n+ay 

If (C2) holds, then 



a" 



1 /■ 1 /, u , 277(a) r 1 , „77(a) 

a J|r|<Q (rl + ay a Jq a^ 

- I 7rT^vi\r\)dr<'- r^dr<2D.^-4. 



□ 



Lemma 6.8 Assume that 77 satisfies (C). Suppose in addition that rj satisfies 
either (CI) or (C2). Then, for any t > and £, S M.'^, we have: 

Nt{0>D['^^=l== [ — -L==-,.{dr), 



(2) 

where D\ is a positive constant depending on t. 

Proof: We denote a = ^^'(^) and b = ^/^TW^■ 

(i) Assume at < 1. We use the fact that cos is decreasing on [0,-|]. Let 
0<c<d<^be such that cos c < cos 1. 

Since the integrand of (j53p is non-negative, the integral is bounded below 
by the integral over the region {c < jrij < d}. For any r in this region, 
(r2 - a2)2 < 2(r4 + a*) < 2{d^ + l)/t^ and 



5? 



{-\ = [cos(ri)-cos(at)]^ > ^ cos^(ai)-cos^(rt) > i cos^ 1-cos^ c M > 0. 
\a/ 2 2 
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Hence, 

NtiO > I rr^-^gl (-) Kdr) > M f f v{dT) A{t). 

■I c<\Tt\<d [T^ - a^r ^aJ ^{t^ + V Jc/t<\T\<d/t 

(68) 



On the other hand, since b> 1, 

Jut' + I ^ Jr (|r| + 1)2-6 (|r| + 6)2 ^ ^ ' 



From ((681) ^^^^ (|69)) . we obtain that: 

('izj ^ssttme at > 1. Using Lemma [6.61 and Lemma [6.71 (with K = l/t), we 

get: 

a2 - Mi/,ay„ (|r|+a)2 ^ ^- Mi/, 6 7r (|t| + 6)2 ^ ^ 

□ 

This concludes the proof of the lower bound and the proof of Theorem 16.21 
The next result is a by-product of the previous lemmas, and gives an alter- 
native condition for the existence of the random field solution to 



Corollary 6.9 Under the conditions of Theorem \6."A for any t > and ^ G M'*, 



1 + *(C) l + vI/(0 ' 

where K^^'' and K^'^^ are positive constants depending on t. Consequently, equa- 
tion ^8\ ) has a random field solution if and only if 



1 + *(0 



^l{d^) < oo. (70) 



Proof: Let a = and 6 = y/l + *(0. The upper bound follows from 

Lemma [631 and Lemma [6.71 fapplied to 6). 

For the lower bound, assume first that at < 1. Then 1 < 6 < \/l + t^. 
Using (|68| and the monotonicity of ry, we obtain that: 

NtiO > A{t) > C{t)^, 

where C{t) is a constant depending on t. 

Assume next that at > \. Then a < b < y/t^+1 a. By Lemma [6.61 and 
using the fact that t] is non-increasing (or 77 is non-decreasing and satisfies (C)), 
we obtain that: 

Nm>Ct^-^>C'{tM, 

6^ 

where C'{t) is a constant depending on t. □ 
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Remark 6.10 In the case of Examples 12.41 and 12.51 with 7 G (— 1, 1), condition 
(|70)) becomes: 

When apphed to Example 1 2 . 4 1 wit h the parametrization 7 = 211 — 1,Hg (0, 1), 
Corollary 16.91 becomes an extension of Theorem 3.1 of W to the case H < 1/2. 
Note that if 7 > (respectively 7 < 0), condition (ITT]) is stronger (respectively 
weaker) than ([551) . 



A The stochastic integral with respect to A4 

Let A4 — {Ai{A); A e TZj,} be a complex random measure on R'' with orthogo- 
nal increments and control measure /I, as specified by Definition 12.21 

Let (p G L^(R'',pt) be a function which is continuous on compact sets. We 
give below the construction of the stochastic integral of (p with respect to Ai . 

Step 1. Jjct A e TZd be fixed. The stochastic integral of (p over with 
respect to M is defined as the L^(Sl)-hmit of Xn = Y^^=iV{'^3,n)M.{Aj^n), 
n > 1, where A„ = {Aj^n)i<j<kn is a partition of A into sets from TZd such 
that ||A„|| = maxjxi;^ l^i.nl ^ Oj Tj^n S ^j,n is arbitrary. The limit exists, 
since the sequence {Xn)n>i is Cauchy: using the orthogonality and additivity 
of Ai, one can prove that E\Xn — Xm\'^ = Iv'n — ^m^^dji, where Lpn = 

^i^j,n)^Aj^,i and |(/?„pdpt — ^ I'jfpc'Mj by the definition of the Stieltjes 
integral (see e.g. p. 228 of 5J. In this calculation, we used the fact that 

E[M{A)M{B)]^Ji{Ar\B), for all A, B e 7^d 

which is again a consequence of the orthogonality and additivity of M. One 
can show that the limit of {Xn)n does not depend on the choice of (A„)„ and 
{Tj,n)j,n- This limit is denoted by MAif) — X4 f{T)M{dT). 
This stochastic integral has^^the following properties: 

(a) Ma{(P + = Ma^^p) + Ma W a.s.; 

(b) EIMaM^ and E[Ma{^)MaW] - Ja vHJl-^ 

(c) M^{ip) = Ma{v), where M-{A) = M{-A) and M{A) ^ M{A). 

Step 2. The stochastic integral of (p with respect to Ai is defined as the 
Lf.(il)-limit of y,i = Ma„{<p), n> 1, where the sequence (A„)„ C TZd is chosen 
such that An C An+i for all n, and U„yl„ = R'^. The sequence (F„)„ is Cauchy, 
since E\Yn-Y,n\'^ = /a,„\a„ I'fil'^djiior any m > n, and /^^ \ip\'^djl^ J^a l^l'^djl 
by the monotone convergence theorem. Since the limit of (Yn)n does not depend 
on the choice of (A„)„, we denote it by A4{ip) = /j^a ip{T)M{dT). This stochastic 
integral enjoys properties similar to (a)-(c) above. In particular, 

E[M{ip)M{tp)] = [ iplpdjl. (72) 
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Moreover, Ai{(p) is a real- valued random variable, for any function <fi which 
satisfies (^(r) = ¥'(— r) for all r e M''. This follows since 

M{ip)= f '^M{dT)= [ ^{-t)M- [dr) = M{ip), (73) 

where we used from property (c) for the second equality above. 

B Some elementary estimates 

Lemma B.l Let 7 € (—1, 1) be arbitrary, (i) For any a > 0, we have: 

I . ^ M -^dT = C^a-^-^, 

where C-y = /^(s^ + l)~^s~'^ds. (ii) For any a> 1, we have: 

C7Wa-^-i< / ^i-^(r2 + l)-7/2rf^<Cf a-^-\ 
Jr T + a 

where C^^ and C^'^ are some positive constants depending on 7. 

Proof: (i) The conclusion follows by the change of variable T/a = T' . 

(ii) We denote by Ii and I2 the integrals over the regions |t| < a, respectively 
|t| > a. We use the notation /(r) ~ ^(t) if c\g{T) < f{T) < c^gir) for any 
T e K, for some constants ci, C2 > 0. When |r| < a, 

(t2 + a2)-i(r2 + 1)^^/2 ^ a-2(r2 + ly/^ ^ a-^r + 1)"^ 
and hence h is bounded below and above by some constants multiplied by: 

/"(r + l)-^dT = a-2-^[(a + 1)^-^ - 1] ~ a'^ia + 1)^"^ ~ a'^'K 
Jo 1-7 

When |t| > a, 

(r^ + a2)-i(r2 + ly/^ ~ r-2(T2 + l)-^/^ ^ r-2(T + 1)"^ ~ r'^-^ 
and hence 72 is bounded below and above by some constants multiplied by: 

Ja 7+1 

□ 
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C Verification of conditions (hyperbolic case) 

Lemma C.l (i) Let ri{T) = |t|^^ with 7 e (—1,1). Then 77 satisfies (C). If 
7 S (0, 1), then rj satisfies (CI). If ^ G (—1,0), then rj satisfies (C2). 

(ii) Letriir) = (1 + t2)-t/2 with j > -1. Then r] satisfies (C). If j G (0,1), 
then 77 satisfies (CI). If J (z (—1,0), then rj satisfies (C2). 

Proof: (i) (C) is clearly satisfied since ri{XT)/?]{T) = X^'* =: C\ for any r > 0. 
If 7 e (0, 1), (CI) holds since for any a > 0, 

/ rj{T)dT — I T^^dr — ' — ari{a). 

Jo Jo 1-7 1-7 

If 7 e (-1,0), (C2) holds since for any a > 0, 

J a 7+1 7 + 1 ^ ^ 



(ii) We first check (C). If 7 > 0, then the inequality 

< C\ (74) 



7?(Ar) / 1 + t2 



77(t) V1 + ^^t' 

is equivalent to 1 + AV^ > C^'^^'^{1 + t^). We choose Ca = {min(l, A2)}-t/2. 

If 7 < 0, then jll]) is equivalent to 1 + AV^ < €^'^^''{1 + r^). We choose 
Ca = {niax(l, X^)}-'^/^. If 7 e (0, 1) then (CI) holds since for any a>K, 



Jo Jo 1-7 Va 



''/^^ a f Ck ^''^ 



1 - 7 V 1 + , 

where = 1 + K^'^. If 7 € (-1,0) then (C2) holds since for any a> K, 



-2 



/•CO 1 

(l + r2)-^/2dr < / r-\CKT^)-^^' ^ Cj.-'/'^a-^-' 

J a 7+1 



^ 7+1 Va / ~ Vl + a 



□ 
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